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Abstract 

We investigate the production of ultra-high-energy cosmic ray (UHECR) in jets from low-luminosity 
active galactic nuclei (LLAGN). We propose a model for the UHECR contribution from the spin- 
down power of black holes (BHs) in LLAGN, which present a jet power Pj < 10'*^ erg s ^ This 
is in contrast to the opinion that only high-luminosity AGN can accelerate particles to energies 
> 50 EeV. We rewrite the equations which describe the synchrotron self-absorbed emission of a 
non-thermal particle distribution to obtain the observed radio flux density from flat-spectrum core 
sources and its relationship to the jet power. In general, the jet power provides the UHECR lu- 
minosity and so, its relationship to the observed radio flux density. We found that the UHECR 
luminosity is dependent on the observed radio flux density, the distance to the AGN, and the BH 
mass, where the particle acceleration regions can be sustained by the magnetic energy extraction 
from spinning BHs and where the strength of the magnetic field at the sites of particle acceleration 
scales with the maximum value of the BH magnetic field, which is ~ 10"^ gauss for a BH of lO^M©. 
We apply the model to M87 and Cen A, two possible sources of UHECRs, whose jet parame- 
ters can be inferred from observational data. Next, we use a complete sample of 29 steep spectrum 
radio sources with a total flux density greater than 0.5 Jy at 5 GHz to make predictions for the max- 
imum particle energy, luminosity, and flux of the UHECRs from nearby AGN. Using our proposed 
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model, it is possible to show that LLAGN can be sites of the origin of UHECRs. In additional, 
the scenario in which the contribution to the UHECR flux from many weak radio galaxies would 
dominate over that from a few strong radio galaxies, or vice- versa, should be substantiated with 
further statistics. 

1 Introduction 



Cosmic rays (CRs) are a direct sample of matter from outside the solar system, and their study can, 
for instance, provide important information on the chemical evolution of the universe or improve 
constraints on Galactic and extragalactic magnetic fields. They can be measured indirectly through 
the study of extensive air showers that are induced as the CRs hit the top of the atmosphere (known 
as CR events). The extensive air showers are currently observed using air fluorescence [e.g.. High 
Resolution Fly's Eye (HiRes) experimen|l| or large array, ground-based detectors [e.g., Akeno 
Giant Air Shower Array (AGASA)i, or both [e.g., Pierre Auger Observatory (Auger)p]. In the 
future, space-based detectors might be another option. UHECR particles are mostly protons or 
fully ionized nuclei with energy above 50 EeV (1 EeV = 10'^ eV). At such high energies, the 
flux of UHECRs is very low and only a few dozen particles per square kilometer per century are 
expected. This is one of the main reasons for the difficulty posed in understanding the origin 
and nature of the UHECRs. Therefore, very large detector arrays are required. The Pierre Auger 
Observatory, by far the biggest cosmic ray detection instrument, uses air fluorescence and water 
detection in a hybrid instrument with an aperture of 7000 km^ sr. 

Joint efforts have been made during the past decade by worldwide, cosmic ray experiments to 
help us understand from where the UHECRs come and what is their nature. It is believed that 
the UHECRs originate in extragalactic sources, as the gyroradius of a proton with an energy of 
100 EeV is of the order of the dimension of our galax y, whereas most of the CR particles w ith 
energy below 50 EeV originate within our galaxy (e.g., Berezinsky et al. 2006; Stanev 2010al lbl). 
If the UHECR particles are protons, they are subject to energy loss by creating pions through 
their occasional collisions with the cosmic microwave background (CMB) photons. This process 
produces a suppression of the cosmic ray energy spectrum beyond 50 EeV, which is known as the 
Greisen-Zatsepin-Guzmin (GZK) cutoff ( Greisen 1966 : Zatsepin & Kuzminlll966 ). Therefore, the 
UHECRs would not be able to survive the propagation from their acceleration sites to us unless 
their sources are located within ~ 100 Mpc. The presence of the GZK cutoff at the expected energy 
in the data released by th e HiRes collaboration was t aken as strong evidence that the UHECR flux 
is dominated by protons (jHiRes -Collaborationll20 1 Oh . 

A suppres sion of the CR flux has also been ob served in the data released by the Pierre Auger 
collaboration (|Auger-Collaborationll2008ai 1201 Oa|) . With respect to primary composition, this col- 
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laboration has exploited the observation of the longitudinal shower development with fluorescence 
detectors to measure the depth of the maximum of the shower evolution, X^ax, which is sensitive to 
the primary mass. A gradual increase of the average mass of cosmic rays with energy up to 59 EeV 
is deduced when comparing the absolute values of Xmax and RMS(Xmax) to air shower simulations 
(| Auger-CoUaborationlbo 1 Oah . 

The present data collected by the Auger Collaboration, which consists o f 69 events of energy 
E > 55 EeV , shows an anisotropy in the arrival direction of the UHECRs (| Au ger-CoUaboration 



2007L l2008bl 1201 Ob|) . Moreover, the arrival direction of the UHECRs is statistically correlated 



with the distribution of nearby extragalactic objects [AGN and gamma-ray bursts (GRBs)], where 
the region around the position of the radiogalaxy Cen A has the largest excess of arrival directions 
relative to the isotropic expectations. 

At highest energies, proton propagation is aff ected only by the CMB, w hereas heavy nuclei 
may be deflected by Galactic magnetic fields (e.g. jMedina Tanco et al.lll998h . 

UHECRs are most probab ly accelerated at astrophysic al shocks, for instance, through a first- 



order Fermi mechanism (e.g., iGallant & Achterberg 



19991) . in very powerful systems that can be 



associated with jets and hot spots in AGN and GRBs. The magnetic field plays an important role 
for the particle acceleration mechanism. The field should be strong enough to confine the particles 
in the acceleration region, but at the same time, weak enough to avoid too much loss by radiative 
cooling. Such shocks can be associated with Poynting flux models for the ori gin of jets f rom fo rce- 
free magnetosphe re above thin accretion disks, which were first proposed by Lovelacel (|l976|) and 
BlandfordI (|l976f) . In the model by Lovelace, the accreting protons are accelerated in the potential 
drop across the accretion disk by electric forces, which then form two high-current, aligned, and 
opposite proton beams. The output electrical power in the beams is proportional to the maximum 
energy of the protons squared, L ~ ^^ax- The maximum energy to which the accretion disk can 
accelerate the proton beams is set by the Eddington luminosity. If one takes the Poynting flux as 
a lower lir nit to the energy flux alon g a relativistic jet, UHECR production in LLAGN cannot be 
explained. iBiermann & et al.l (|2008|) rewrote Lovelace's equation as 



- —/flare , „ 



(1) 



where Z is the mass number of the nuclei, y^^ is the Lorentz factor of the shock, and /flare(< 1) is 
the intermittency. As the authors state, probably all three elements are required if one considers 
UHECR production by sources like M 87 and Cen A, whose energy fl ow along the jet are < 10"^^ 



erg s ' and < 10"^^ erg s \ respectively (IWhysong & Antonuccill2003|) . 



As an alternative, iFarrar & Gruzinovl (|2009|) showed that very intense, short-duration AGN 
flares that result from the tidal disruption of a star or from a disk instability can accelerate UHE- 
CRs. On the other hand, magnetic reconn ection in relativistic jets represents another option for 
UHECR acceleration (e.g.. lGianniosll20 1 01) . 
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Boldt & GhoshI (119991) suggested that particles with energies > 10 eV may be accelerated 



near the event horizons of spinning BHs associated with presently inactive quasar remnants. The 
required electromagnetic force is generated by the BH induced rotation of externally supplied 
magnetic field lines threading the horizon, where the magnetic field is supported by external current 
and the horizon is an imperfect conductor with resistance ~ 100 f2. Therefore, the BH behaves as 
a battery, driving currents around a circuit, with an electromagnetic force of up to 10^' eV for a 
BH with a mass of 10^ Mq (e.g., lLovelacelll976l) . In this case, the production of observed flux 
of the highest energy cosmic rays would constitute a negligible drain on the BH dynamo. That 
is, replenishing the particle ejected at high energies (> 10^° eV) would require a minimal mass 
input; a luminosity of lO'*^ erg s^' in such particles (if protons) corresponds to a rest mass loss 
< IO^^Mq in a H ubble time. Particle acceleration to UHEs from the spin-down power of BHs was 
also discussed in Blandford where the acceleration regions can be sustained by magnetic 

energy extraction from spinning BHs. 

In this paper, we propose a model for UHECR contribution from the spin-down power of BH 
in LLAGN. The particles in the jet manage to tap the spin-down power of the BH and then are 
accelerated at relativistic shocks with energies up to the UHE domain. The electrons lose their 
energy through synchrotron emission, whereas the protons are capable of surviving the radiative 
cooling and perhaps of propagating through the intergalactic and Galactic medium towards us. 
Since both particles undergo the same acceleration process, there must be a correlation between 
the electron synchrotron emission and the UHECR proton energy. We seek this correlation to make 
predictions for maximum energy, luminosity, and flux of the UHECRs from nearby LLAGN. This 
is in contrast to the op inion that only high-luminosity AGN can accelerate particles to UHE domain 
(e.g.. lZawetaDl2009l) . 

In Section [2l we provide a description of the model. In Section [3l we derive the relation 
between the jet power and the observed radio flux density for a flat-spectrum core source. Based 
on this relation, we derive in Section H] the luminosity and flux of the UHECRs. In Section [51 
we calculate the particle maximum energy taking into account the spatial limit and synchrotron 
emission losses. In Section [6l we present the application of the model to M87 and Cen A. Both 
sources are LLAGN with a mass accretion rate relative to the Eddington accretion rate less than 
~ 10"^ (therefore, they can be powered by the BH spin down) and present strong radio-emitting 
jets. In Section|71 we provide the predictions for nearby galaxies as possible sources of UHECRs. 
In Section [81 we present a summary of the key points and discuss the implication of this model for 
further studies of the UHECRs. 
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2 Model description 



2.1 Model conditions 

• We assume that the UHECRs are accelerated by shocks in AGN jets, which are launched 



from the inner accretion disk which is located inside the BH ergosphere (IDutan & Biermann 



20051) . The inner disk extends from the stationary limit rgi inward to the innermost stable 
orbit Tms- When the mass accretion rate relative to the Eddington accretion rate is ~ 10"^, 
the jets can be powered by the spin down of the BH, which is transferred to the disk inside 
the ergosphere by closed magnetic field lines that connect the BH to the accretion disk. The 
jet propagates along a cylinder of length zo (see Fig. [B using the BH spin-down power and 
then extends into a conical shape with a constant opening angle 2 0, as a consequence of the 
free adiabatic expansion of the jet plasma. (T he tip of the cone is lo cated at some z < zo.) A 



similar geometry of the jet was considered by iMarkoff et al.l (|2001|) 



The calculations are performed for the case when the UHECRs would have been protons. 
For heavy nuclei of a given atomic number (Z), the particle energy will scale up with Z. 

In the observer frame, the magnetic field along the jet varies as 5 ~ T, ^i^d the electron 
number density in the jet scales as ~ 7jZ"^, where yj is the bulk Lorentz factor of the jet. (See 
discussions below.) 

We set the slope of the particle density distribution to p = 2, which corresponds to a flat- 
spectrum core sources with a spectral index a = 0.5, and the strength of the BH magnetic 
field to its maximum value (B^^). The latter condition provides, in turn, the minimum 
values of the particle maximum energy, luminosity, and flux of the UHECRs. In Section [6l 
we discuss some implication of a different choice for the value of p. 



2.2 Magnetic field scaling along a steady jet 



To describe the jet physics, we use the following reference frames: (i) the frame comoving with 
the jet and (ii) the (rest) frame of the observer, in which the relativistic jet moves with the bulk 
Lorentz factor. 

In a frame comoving with the jet, the poloidal component of the magnetic field is considered 
to vary as Bp ~ z~^. This variation follows from the conservation of magnetic flux along the axis z. 
To keep the field divergence -free, the toroi dal component must vary as 5^ ~ z'^. This topology of 



B^ ~ z ^ was first derived by iParkeii (|1958|) for the magnetohydrodynamics solution of a spherical- 



symmetric flow (so that, a j et can be considered a conical cut along the flow surfaces). [See also 
Blandford & Konigll (119791) . 1 The observational support to this variation of B^ is specified later 
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Figure 1 : Schematic representation of the jet geometry. The jet is launched from the inner disk, 
extending from the stationary limit inward to the innermost stable orbit, and propagates along a 
cylinder up to a distance of zo- It then expands freely into a conical geometry with a constant 
opening angle 2 6. The magnetic field lines threading the disk near the BH (dashed lines) are 
wound up, far from the BH, into a toroidal magnetic field that coUimates the jet. 



on in this section. At a distance, say, zo, the poloidal and toroidal components of the comoving 
magnetic field become approximately equal 5po - fi^o- We consider zo of a few gravitational radiijj 
based on the fact that the VLBI observation, for instance, of the jet in M87 at 7 mm gives evid ence 
on the jet coUimation (by the toroidal magnetic field) on scales of 60-200 (|Biretta et al. Il2002h and 
the gl obal 3.5 mrn VLBI observations have resolved sizes for the compact radio sources of ~ 10 
(e.g., iLee & et al.ll2008|) . A large-scale and predominantly toro idal magnetic field can exert an 
inward force (hoop stress), confi ning and coUimating the jet (e.g. jBisnovatyi-Kogan & Ruzmaikin 
19761: iBlandford & Paynelll982|) . The magnetic hoop stress is balanced either by the gas pressure 
of the jet or by centrifugal force if the jet is spinning. From zo upward, the poloidal component 
of the magnetic field becomes weaker, so that the field lines are soon wound up in the azimuthal 
direction by the jet rotation. Thus, above zo, the magnetic field along the jet is nearly azimuthal 
B ^ Bip (for a steady jet) and varies inversely proportional to the distance along the jet: 



B = Bo\- 



(2) 



where Bq = B^q ^ B^q is the strength of the magnetic field at the height z = Zq above the equa- 
torial plane of the BH. This z- dependence of the ma gnetic field appears to be contradicted by 
radio-polarization observations ([Bridle & Perleylll984|) . These observations strongly suggest that 



^The gravitational radius is defined as = GM/c^ = r^(M/10'^Mo) = 1.48 x lO'^^CM/lO'^Mo) cm, where G is the 
Newtonian gravitational constant, M is the BH mass, and c is the speed of light. 
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the magnetic field is predominantly parallel to the jet axis initially and only later becornes per- 
pendicular to the jet axis, with some parallel magnetic field left over. iBecker & BiermannI (|2009h 
argued that the basic pattern of the magnetic field is indeed 5^ ~ z ' and that the observational 
evidence for a parallel magnetic field is due to highly oblique shocks. Their argument is based on 
the observations of the jet structure which indicate the occurrence of the moving shocks between 
20 and 200 r^, w hile the first stationary, strong shock can be produc ed in the approximate range of 



(3 - 6) X 10^ rg (iMarkofFet al.ll200lL 120051: iMarscher & et al.ll2008h 



The strength of the magnetic field i n the comoying f rame Bq can be related to the poloidal 
magnetic field in the BH frame Bu (e.g.. lDrenkhahnll2002r) as 



H 



1 5;:;'"' / b 

Bn = — B\i — I 



where the maximum value of the BH magnetic field is given by 

,-1/2 



(3) 



5^?^" ^0.56x 10'(— ^ 



gauss. 



(4) 



which is obtained in a similar manner as the calculation performed by IZnajekl (|1978|) . with the 
diff"erence that we set the BH potent ial drop to the specific energy of the particles at the innermost 
stable orbit, whereas IZnajekl (|l978|) makes use of the fact that the Eddington luminosity sets an 
upper bound on the radiation pressure (as the disk is radiatively efficient). The maximum value of 
the BH magnetic field corresponds to the time when the accretion rate was as high a s the Eddingto n 
accretion rate. In this case, the BH spin parametej^ is limited to a* = 0.9982 ( Thome 1974 ). 
Although this limit might be even closer to the maximal value of the spin parameter ~ 1, this will 
introduce just a small variation of the maximum value of the BH magnetic field. 



2.3 Electron and proton number densities 

The jet is assumed to be composed mainly of electrons, positrons, and protons. We denote by 
/ep = He/Hp the ratio of the electron to proton number densities, where the number densities are 
measured in a frame comoving with the jet plasma. Unless otherwise noted, should be assumed 
to include the positron number density as well. It is straightforward to generalize to a mixed 
chemical composition, including many heavy nuclei. Furthermore, both electrons and protons 
can have thermal and non-thermal populations before being accelerated at the shock. There may 
also be a substantial number of positrons from pion production and decay processes (also called 
secondaries). 

*The BH spin parameter is defined as a, = J/Jma\(- a/r„), where a - J/Mc is the angular momentum of the BH 
about the spinning axis per unit mass and per speed of Ught and /max - GM^/c is the maximal angular momentum of 
the BH. Furthermore, the BH spin parameter obeys the condition: -1 < a* < +1. 
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Now, we look for the expression of the proton and electron number densities injected into the 
accelerating region. First, we consider the mass flow rate into the jets, which in the comoving 
frame is given by 

^ dt ^'^'^'^ ^ It ["'^^^'^^-■=0] = «mvj(5)-=o, (5) 

where pj is the rest-mass density of the jet, Vj is the comoving volume of the jet, (5 )-=o is the 
launching area of the jet, z is the length of the cylinder along which the jet propagates before 
expanding freely in a conical geometry, and vj = ySjc is the bulk velocity of the jet. 
The surface area between two equatorial surfaces of a Kerr BH can be calculated as 

(J5),=o = (^) 2ndr, (6) 

where the Kerr metric functions are: 

A = - Ir^r + and A = / + r^a^ + Ir^ra^, (7) 

where r is the coordinate radius. Next, we use normalizations to the gravitational radius, so that 
= r/r„ is the dimensionless radius. The surface area is then: 



I 



where the factor increases from ~ 2 to ~ 80 as the BH spin parameter increases from 0.95 to ~ 1 . 
For the first equality, we use the fact that the inner disk, from where the jet is launched, has its inner 
and outer radii at the innermost stable orbil0 Vs and stationary limit = Ir^ = tq, respectively. 

The comoving density of the jet can be expressed in terms of the ratio of the electron to proton 
number densities: 



m, 



nm = npHip + n^m^ = n^m^ 1 1 + /ep— j = n^m^fo. (9) 



For protons dominating over the electrons, /gp < 2x10^, where electrons and positrons can partially 
occur as secondaries. This fraction of CR is ~ 10"^ from data of CRs at 1 GeV. One can get 
a ratio of unity assuming tha t the spectra go down to rest mass, which is implausible [see, e.g.. 



Protheroe & BiermannI (|1996|) and references therein]. 

Substituting Eqs. ([8]) and ^ for (|5]), we obtain the mass flow rate into the jet in the observer 
frame by including yj: 

Mj = y^fripmpfolTTrlko., (10) 



^Once the accretion flow reaches the innermost stable orbit, it drops out of the disk and falls directly into the BH. 
The expression for the radius of the innermost stable orbit r^s is given by eq. (2.21) in .Bardeea (,1970i) . 
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This expression provides the proton number density, which we use to derive the electron number 
density: 

yjfifmpfolnrjko 

We shall use this result later for evaluating the self-absorbed synchrotron emission of the jets 
(Section 1231). 



2.4 Particle energy distribution 



We suppose that a shock is produced at the jet height z = Zo (e.g. jMarkoff et alj|2001|) . As a result, 
a power-law energy distribution of the particles is established. For a given frequency, the emission 
from the synchro tron process for electrons and protons gives a factor (mp/mp)^ ^ 10'" in favor 



of electrons (e.g.. iNovikov & Thomdll973|) . In addition. iBiermann & Strittmatterl(ll987|) showed 
that the proton synchrotron emission can be competitive if one considers that the proton emission 
ranges to much higher photon energy. Moreover, it is also not at all obvious that they have the same 
normalization at the same Lorentz factor of the particle and that the particles have a continuous 
power-law from rest mass to UHE. The number density of the electrons in the energy interval E, 
E + dE [or rfieC^y, m^c^{y + dy)] has the power-law form: 



N{E)dE = CE-PdE 



(12) 



or, in terms of the Lorentz factor. 



N(j)dy = Cy-^dy, (13) 
where y e rymiii,ymax1 is t he Lorentz factor of the electrons and p is the power-law index (e.g.. 



Rybicki & Lightmanll 19791) . The normalization coefficients of the electron number density in Eqs. 
[T2l and [13] are related by 



C = C'im^cy 



(14) 



The normalization of the electron number density C follows the pressure in a reheating flow 
(~ z~^), where as an adiabatic flow would give a steeper dependence, which leads to shocks (e.g.. 



Sanderslll983b . Adiabatic behavior implies P ~ p^"', where y^^ is the adiabatic index. Since a rela- 



tivistic fluid usually has yad = 4/3, in an adiabatic flow the temperature runs as Furthermore, 
one can have a conical flow only if the temperature of the flow is approximately constant. The 
energy for reheating can be taken from the flow through highly oblique shocks. The normalization 
of the electron number density, in the case of a conical jet, is therefore: 



(cm-^). 



(15) 



Possible values of the power-law distribution index p of the electrons accelerated by the rela- 
tivistic shock are discussed later in Section |6] 
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2.5 Self-absorbed synchrotron emission of the jets 



The spectra from compact radio sources can be explained by self-absorbed synchrotron emission of 
the jets produced by electrons with a power-law energy distribution. In this s ection, we rewrite the 



quanti ties which describe the self-absorbed synchrotron emission, derived in lRybicki & Lightman 



(Il979h . and express them under the considerations of the model presented here. We first introduce 



the absorption coefficient, optical depth, synchrotron emissivity, and source function in order to 
calculate the flux density of the synchrotron emission from radio sources with a flat-spectrum 
core. (The quantities which describe the self-absorbed synchrotron emission are in cgs units.) 

In a frame comoving with the jet plasma, the absorption coefficient of the synchrotron radiation 
can be calculated as 

_ V3e^5 sinao /3e5 sinQ'o\^ ^ 

8mne \ 2mnlc^ ) ^^^^ 

3d -I- 2\ ^ /3p -I- 22\ ^;-+4 

' • ^ (cm ), 



12 / \ 12 

where e is the electron electric charge, is the electron mass, p is the power-law index of the 
particles distribution, C is the normalization factor for the power-law electron energy distribution 
(Eq. [I2l), B is the magnetic ffeld in the frame comoving with the jet, Y(x) is th e Gamma function of 



argurn ent x, and y is the frequency of the synchrotron radiation [eq. 6.26 in Rybicki & Lightman 



(|1979[) 1. The average over the pitch angle a^, for a local random ly oriented mag netic field with 



a probability distribution ^ sin q'o'^^Q'o, is given by the integral in iLongairl (|1994l) . Including the 



values of the physical constants and using the expressions for the normalization of the electron 
distribution function, C and C (Eqs. [14] and [15]), as well as for comoving magnetic field along the 
jet (Eq.l2]), the absorption coefficient then becomes: 

ay = K,CJ-\ B^-v--, (17) 



where 



= 8.4 X 10~-^(1.25 X 10''^)5 (8.2 x 10"^)^' ' — 



(18) 



To calculate the observed distance along the jet where the jet becomes self-absorbed, we first 
determine the optical depth Tv of the jet material. The averaged path of a photon through the jet has 
the length r(z), which is a reasonable approximation for a jet observed at large inclination angle 
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(e.g.. lKaisej|2006f) . We introduce a factor Iq in the expression of the path length to account for a 
small inclination angle. Thus, we can write the optical depth as 



Ty = ayr{z)lo. 



(19) 



For conical jets, the intrinsic half-opening angle is given by tan^ = r/z = ro/zo- With the absorp- 
tion coefficient specified through Eq. (fTTT ). the optical depth can be written as 



Ty = KiCf^rolol — ] ^ 



j,+4 



Ell _p+4 



(20) 



One can define the distance along the jet where the jet becomes self-absorbed Zssa as the dis- 
tance z for which Ty = 1. Using Eq. (I20l) . one obtains: 



p+6 £+2-1 
1 /;+4 jy p+4 



(21) 



The total power radiated per unit volume per unit frequency by a non-thermal particle distribu- 
tion equals: 



smao / meCco 



InrUf-c^ p + l \3e5sinao 



2 



19 



1 



(22) 



where oj = 2nv [Eq. 6.36 in lRybicki & LightmanI (Il979l) 1. Using Eq s. Q and (fT5l) , as well as the 
method to calculate the averaged pitch angle employed in iLongaii the expression of the 

total power becomes: 



Py = 2nP„. = K2C'A- 



£+5 

2 £+1 ,,-1 



(23) 



where 



ii:2 =3.7x 10-2^(l.2x 10*^)"'' (p + 1)"'^ 



^'?-nn?-T2ir 



p -I- 5 
4 



— 1 



P + 7 
4 



(24) 



Next, we derive (z/zo) from Eq. |20]when Ty = 1. With this, the expression for the total power takes 
the form: 

(25) 
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and the emission coefficient is simply jV = Py/4-n. 

The emission coefficient is defined as the product between the absorption coefficient ay and 
the source function Sy. At low frequencies, the emitting region is opaque to synchrotron radiation 
and the observed intensity of radiation ly is proportional to the source function, while at high 
frequencies, the region is transparent and the observed intensity is proportional to the emission 
coefficient. The two dependences should be matched at the transition from opaque to transparent 
regimes. This transition corresponds to an optical depth Ty = 1. The source function in the self- 
absorbed limit is then: 



S. = ^- = ^^ {-] B-K^ (erg s-'cm-^Hz-i), (26) 
An ay AnKi \zol 

where the last equality was obtained using Eqs. (fTTl) and (1231 ). Note that in the source function, 
the emitting frequency does not depend on the power-law index of the electron energy distribution. 
For Ty = 1, the source function becomes: 

S; = K, (C^ro/o)^ Bl^'v" (erg s-^cm-^Hz-'), (27) 

To obtain the emission spectrum, one needs to solve the equation for the radiative transfer 
through a homogeneous medium. Because the angular sizes of the jets are small, instead of the 
specific intensity of the radiation, one usually measures the flux density Fy (energy per unit time, 
per unit frequency interval, that passes through a surface of unit area). Thus, 

dFy = lydQ. = S y [ 1 - exp(-Ty)] dQ. (28) 

Because the frequency shift of the approaching photons, specified by the Doppler factoi@, is Vobs = 
DjV, the transformation of the specific intensity to the observer frame is: 



' v,obs 



n'ly, (29) 



where the relativistic invariant quantity ly/v^ was used. The solid angle corresponding to the source 

dn = (30) 



where is the distance from the observer to the jet source and r = z tan 6. If we insert Eq. (|30l) 
into Eq. (l28l) and integrate it from zq to z, we obtain the flux density of the synchrotron emission 



'^The Doppler factor of the jet is £)j = - /Sj cos^)"', where ip is the inchnation angle of the jet axis with 

respect to the line of sight (which is Lorentz transformed through sini^obs = sin 1/3). The angles are rotated by the 
Lorentz transformation, so that a jet seen at angle yr' is rotated basically to a transverse view for large yj. 
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in the case of = 1 as 



F; = 5;[1 - exp(-l)];r(tane)D-V 



(31) 



where the second term in the last squared bracket can be neglected with respect to the first term for 
Zo (where the jet emission becomes self-absorbed). Using Eqs. (|2TI) and (ITT]) , the flux density 
is then: 

^ 2p+13 2p+3 

(32) 



^ 2p+13 2p+3 

F' = i^4(QVo)^V^' B^'-' D;\tanerK 



where K4 = 0.16^^ ''^'^K2. The radio flux density in Eq. (|32l ) does not depend on the emitted 
frequency of the radiation since we already adopted the case of flat-spectrum core sources when 

Tv= 1. 

For a power-law synchrotron spectrum of the form Fobs 
related to the intrinsic flux density as 



Fobs = £>f"F', 



v^^^, the observed flux density is 

(33) 



where Eq. (1291 ) was used. Equation (1371 ) is valid for a single blob emission (e.g. jBlandford & Konigl 



19791) . For a continuous jet consisting of uniformly-spaced blobs, one needs to consider the emis- 
sion per unit length in the observer frame. Therefore, 



Fobs = "F', 



(34) 



due to the fact that the number of the blobs observed per unit length is ~ 1 /£)j. 



3 Relation between the jet power and the observed radio flux 
density for a flat-spectrum core source 

In the previous section, we established the expression for the radio flux density from flat-spectrum 
core sources (Eq. |32]). This quantity reflects the r adiative property of the jet, as the radiated energy 
is replaced by dissipation of the jet power (e.g., iBlandford & Konig]||l979l) . In this section, we 
seek the relation between the jet power and the observed radio flux density. First, we consider the 
jet power in the observer frame defined as 



(35) 



which follows, e.g., from iFalcke & Biermann (Il995b [see also IVila & Romerol (|2010l) 1. and for 
which we need to evaluate Mj using Eq. \TT\ An upper limit for the electron density is specified 
by He < C(). So, we can substitute Eq. (fTTI) for the expression of the observed radio flux density 



13 



(either Eq. |33]or Eq. |34l depending on the structure of the jet) and find the mass flow rate into the 
jet Mj. The strength of the magnetic field Bq follows from Eqs. Q and dH). This procedure yields 
the power of the jet: 



2;)+ 13 ^ Pj± 



0.05 



2r„ 



2p+13 
S 



Dmax 



where 



n 



2p+3 

5 £+4 2(/)+4) 

^obs 



.Ell 2p+13 



M 



lO^Mp 



2p+3 
10 



(36) 



ergs 



-1 



K5 = -m^c'K^ ' (5)^(2.96 x lO'^)"^ 
(0.56X 10^^(0.05)"^^, 



p+13 



(37) 



where h = {{p + 5){p + 4)]/10 for single blob emission (Eq. \^ ov h = {{p + 3)(p + 4)]/10 for 
continuous blob emission (Eq. [341) and / = fokoQofepT^ ■ We use a normalization value for the 
Lorentz factor of the jet, say 5, although this factor can range from ~ 2 to ~ 100, as observational 
data suggest. We adopt /gp ~ 10"^ (and then /o - 1), which means that there is, in average, 
one hundred electrons/positrons for at least one proton in a jet that is powered by a very rapidly 
spinning BH (a* > 0.95) and ob served at a large angle (> 1 0°). Flat spectrum cores are predicted 
for any angle to the line of sight (|Blandford & Kdnigllll979t) . They are pointing close to the line of 
sight only if the cores dominate over the extended emission. 



Mildly-relativistic shocks (1 < 7s ^ 30) are believed to occur in the AGN jets. iGallant & Achterberg 



(| 19991) have shown that at relativistic shocks, the particles typically perform a fraction ~ y^, of 



the Larmor orbit upstream before recrossing the shock, and the particle energy gain increases by a 
factor of y^jj in the first shock-crossing cycle. In the subsequent shock-crossing cycles, the energy 
gain is of the order of 2. The predicted spectral index is =^ 2. A series of Monte Carlo simula- 
tions (which are meant to find a way of constructing the trajectories of particles whose distribution 
obey the desired transport equation), performed under a wide range of background conditions at 
the shock front, indicate a value p - 2.2 - 2.3 for the slope (e.g ., iBednarz & Ostrowskilll998 : 
Achterberg et al.lboOli iKirk et al.lboOOl: iKeshet & Waxmanll2005|) . It follows that in the regime 
of arbitrarily high Lorentz factor shocks, the acceleration process generates particle spectra which 
are quasi-independent of the considered background conditions, leading to a quasi-universal slope 
of ~ 2.2. This picture can slightly be changed when one considers more realistic conditions in 
the vicinity of the shock. For example, iNiemiec & Ostrowskil (|2006h studied possible models for 
perturbed magnetic field upstream of the shock and found that for superluminal, mildly-relativistic 
shocks, a flattening of the spectrum occurs, p - 1.5, with a cut ofl" at lower energies than as 
expected for UHECRs. A flattening of the spectrum, p ~ 2.1 - 1.5, is also observed in simu- 
lation results obtained by iMeli et al.l (|2008|) for the case of superluminal, ultra-relativistic shocks 
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(100 < ysh < 1000), with a turbulent magnetic field and various shock obliquity; however, for 
subluminal, mildly-relativistic shocks, the spectral slope has values between 2.0 < p < 2.3. 

For the following calculations, we adopt p = 2, which is the upper value of the spectral index 
for flat-spectrum core sources. The jet power (Eq. |36j then becomes: 



Pj =8.78 X \0^%DJ^ (y) 

yioaJy) \Mpc 



Dmax 



tan^y 
005/ 

Apcj \ 



6/5 



2r„ 



M 



lO^Mp 



-17/5 



-7/10 



(38) 



where a continuous jet was considered (h = 3). For p < 2, the results for the jet power are slightly 
reduced. 



Blandford & Konigll (119791) showed that the flat-spectrum radio synchrotron emission of a com- 



pact jet core is produced by superposition of self-absorbed synchrotron spectra at diff'erent posi- 
tions of the jet. In their model, the observed radio flux density depends on the jet power and the 
distance to the jet source as 

f,^,-p"''^d;\ (39) 



Heinz & SunyaevI (120031) obtained a generalization of Eq. (|39l ). for any scale-invariant jet 



pin+Sa)/l2 j^- 



model producing a power-law synchrotron spectrum with an index a in the form of Fobs 
where M is the BH mass. Because of the large mass difference between AGN and microquasars, 
this non-linearity function of the radio flux density with the BH mass indicates that the AGN jets 
are more radio-loud tha t the rn icroquasar jets. 

Falcke & Biermann (1995) found that for radio-loud|^ AGN, the observed radio flux density 



depends non-linearly on the BH mass: 



obs 



^1.42^2.2^-1. 



(40) 



This result is obtained for the case of an accretio n-dominated jet, Pj ~ Ldisk, where Ldisk is the 
luminosity of the disk (see also lFalcke et al.lll995|) . 

Using Eq. (|38] ) for the observed radio flux density of a conical jet, we obtain: 



Fobs - P]"'D:'-M'l''D]'^y:''l\i&ner\ 



(41) 



where the magnetic field along the jet varies as 5 ~ yj z , the electron number density in the jet 
scales as ~ yjz"^, 5h - 5™^", and = 2 r^. Since the observed radio flux density in Eq. (1411 is not 
dependent on the distance along the jet, the expression can be applied to microquasars as well. 



'The radio-loudness parameter is defined as the ratio of radio luminosity emitted by the jet to the UV/X-ray 
luminosity emitted by the accretion disk. 
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4 Luminosity and flux of the UHECR 



In this section, we seek for the UHECR luminosity flux (Fcr) specified as a function of the ob- 
served radio flux density. First, we consider the UHECR luminosity defined as 

^CR = ecRPj = ecRTiMf^, (42) 

where it is assumed that the UHECR luminosity is a fraction (eca) of the jet power, with Pj = 
^kin "I" ^magn LcR. If we were to adopt the point of view that the jet power is shared equally 
in a comoving frame between the baryonic matter, magn etic field, and cosmi c rays extending to 
the highest energy, £cr - 1/3. In the jet-disk model of iFalcke & BiermannI (| 1995b . the energy 
equipartition in the comoving frame appears to be a good approximation. It would also suggest 
that AGN driven by the BH spin-down power, and sufi'ering from a low mass accretion rate, may 
attain a higher Lorentz factor, consistent with some observations. 

Using the expression for the jet power (Eq. |38), the cosmic ray luminosity (Eq. I42l) becomes: 

17/5 /tan^x6/5/ ^0 



X' (D^r i 



\mJy/ \Mpc/ \\QPM^ 
Given the UHECR luminosity, we can easily obtain the UHECR flux: 

^CR 



CR 



AnDV 



(44) 



where we do not include the cosmological distance as we refer to nearby radio, flat-spectrum core 
sources with a redshift up to z ~ 0.018. 



5 Maximum particle energy of the UHECR 

Now, we look for the maximum energy of the UHECR in the case of the spatial (geometrical) limit 



(IFalcke & Biermannlll995i): i .e.. the jet particle orbits must fit into the Larmor radius. Conform to 



Gallant & AchterbergI (|1999|) . the maximum particle energy in the downstream rest frame can be 



written as 

£max = YseZBor, (45) 

where 7s is the Lorentz factor of the shock and Z is the particle mass number. Using the expression 
for the magnetic field along the jet (Eqs. |2]and|3]) and the fact that tan 9 = ro/zo, the maximum 
energy of the UHECR particles (in the observer frame) becomes: 

EZ. = eZBrro[^^[-^). (46) 



16 



For protons, 

1/2 



^X. = 5x.0^»(|l)(|l)(_^) ,eV). (47) 

where js - yj was used. 

Next, we look for the maxim um energy of the UHECR in the case of the synchrotron loss limit 



(IBierrnann & Strittmatteiill987j). Sett ing synchrotron losses equal to diffusive shock acceleration 



gains. IB iermann & Strittmatterl(|1987l) showed that a ubiquitous cutoff in the non-thermal emission 
spectra of AGN can be explained. This requires that the protons to be accelerated near 10^' eV. 
The frequency cutoff (Vh.) might be pro duced at about (3 - 6) x 10^ rg. R ewriting the expression for 



the maximal proton energy derived by iBiermann & Strittmatteii (|1987l) . 



E^Zx ^ 1.4 X 10^° ( \ — B-^l^ (eV), (48) 

l3xl0^4Hz/ ^ ^' ^ ' 

and using the expression for the magnetic field along the jet (Eqs. [2land|3]), the maximal proton 
energy in the loss limit reads: 

Is X 1014Hz/ 15/ 

D \-l''2 / ^1/2/ vl/2 (49) 



6 Application to M87 and Cen A 

In this section, we investigate the UHECR luminosity flux for two possible sources of UHECR, 
M87 and Cen A, whose jet parameters can be inferred from observational data. 

Taking r = Ir^ and Bu ^ 5™™, the equations for the maximum particle energy in the spatial 
(Eq. 1471) and loss (Eq. |49l) limits, as well as for the UHECR luminosity (Eq. I43l) . become: 

^ 4.2 X 10'^ ' * ^ 



3 X 1014Hz/ I 5 



- (eV), (50) 



= 3.65xl03^,(l-,cos.f(|r(^f 
mJyj \Mpc/ \10''Mg 



.6/5 j vl2/5 I v-7/10 

(ergs ). 
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This set of equations is used for the following estimations on UHECR contribution from the BH 
spin-down power. For the expr ession of E^^^^ in Eg . (l50l). we use (z/zo) ~ 10^. Our choice is 
based on the results obtained by iBecker & BiermannI (|2009[) . which show that a first large stead y 
shock can be produced at about z ~ 3 x 10^ rg [following the work by Markoff et alJ (|200lh l. 
This is confirmed by ob servations of a blazar inner jet as revealed by a radio-to-y-ray outburst 
( Marscher & et al.ll2008|) . The same conclusion was reached earlier by iBiermann & Strittmatter 
(|l987|) using the observed cutoff" in the radio emission of AGN. 



Table [U contains our estimations for the maximum particle energy, luminosity, and flux of the 
UHECRs in the case of M87 and Cen A, whose jet parameters can be obtained from observational 
data. The observed radio flux density of the core corres ponds to a frequency o f 5 GH z. For 
comparison, we use the energies along the jet estimated by IWhysong & Antonuccil (|2003h . which 
are ~ 10"^^ erg and ~ 10"^^ erg for M87 and Cen A, respectively. The estimation of the 
luminosity and flux of the UHECR corresponds, however, to the upper limit of the slope of the 
particle density distribution of p = 2. For a steeper slope of the particle density distribution, the 
luminosity and flux of the UHECR increase, whereas for a flatter slope they decrease. For instance, 
if we take p = 2.4 (so that a = 0.7, the very low limit for a steep synchrotron emission spectrum), 
the luminosity and flux of the UHECR increase to values which are ~ 2.71 times those when p = 2 
for M87 and ~ 3.38 for Cen A. Although both sources have a low jet power, they are powerful 
enough to provide the environment for particles to be accelerated to ultra-high energies of ~ 1 ZeV. 
For both sources, the jet power is supplied by the BH spin-down power, as their mass accretion 
rate relative to the Eddington accretion rate is less than ~ 10"^. 



7 Predictions for nearby galaxies as ultra-high-energy cosmic 
ray sources 



We apply the model pr oposed here to a complete sample of steep spectrum radio sources ([Biermann & et al 



20081 : lCarametell201 11) . at redshift z < 0.025 (about 100 Mpc), with a total radio flux density larger 
than 0.5 Jy. The numbers for the estimated flux and maximal energy exclude the GZK eff"ect, but 
includes the distance efi'ect. The selection criteria used by the authors are presented in more detail 
in their papers. Table |2]lists the predictions for the UHECR particle maximum energy, luminosity, 
and flux. We emphasize that there could be a common scaling limit, such as a condition that the 
Larmor radius has to fit three times or five times into the jet. The scaling limit is not critical to our 
predictions as long as we refer the quantities to, say, those of M87; therefore, the jet parameters 
are assumed to be the same as for M87, all scaled by BH mass, radio power, and distance to the 
source. This can be seen by comparing the estimated values of the UHECR maximum particle 
energy and flux for Cen A in Table [i|(with the jet parameters inferred from observational data) and 
Table[2](with the jet parameters of M87); the differences are within one order of magnitude. Using 
a scaling relative to M87, one implicitly allows for the possibility that all sources produce higher 
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Figure 2: Integral source counts of the same sample in Table |2l The model parameters for the best 
fit of a power-law model (N = C f) with ^(> 0.405) are: C = 4.88 + 1 .74 and k = 0.86 + 0.45. 

nuclei (Z > 1). Of course, this again assumes that all sources are the same in this respect. We ar- 
gue that although the sources are LLAGN, they can be sites of accelerating particles to ultra-high 
energies ~ 1 ZeV with an UHECR luminosity < lO"^"* erg s~K 

In Fig. [21 we show the plot of integral source counts (N) versus relative UHECR flux (last 
column in Table O, where the data points are split in 8 bins of variable width with the bin upper 
edge three times larger than the bin lower edge. Although the data points are too sparse to dis- 
criminate models of the flux distribution, we estimate the most likely model parameters that could 
generate these data points. We fit a power-law model (A'^ = C f) to the mean values of the data 
in each bin for A'^(> 0.405) using the log-likelihood method0 The errors for the data points are 
estimated by VA^. For the best fit that we obtained, the model parameters are: C = 4.88 ± 1.74 
and K = 0.86 + 0.45. (Since the source counts for each bin is too low, especially in the tail of 
the assumed power-law distribution, we cannot perform a ;^f^-test.) When k < I, the contribution 
to the UHECR flux from a few strong radio sources dominates over that from many weak radio 
sources. On the other hand, within the errors we can also consider a power-law model with a: > 1, 
which implies that the contribution to the UHECR flux from many weak radio sources is domi- 
nant. To distinguish one possibility from the other one, more information on radio data and BH 
mass from other nearby AGN, in additional to those listed in Table|2l are required; so that, a larger 
complete sample of nearby AGN may lead to an improved statistical model to answer the question 
of whether the flux of the UHECRs is produced by many weak radio sources or a few strong radio 
sources. 

"^To plot numerical results and to analyze data, we use ROOT, which is an open source, data analysis program 
developed at CERN. 
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8 Summary and conclusions 



In this work we developed a new model for UHECR contribution from BH spin-down power. We 
relate the observed radio flux density to the luminosity and flux of the UHECRs, and calculate the 
maximum particle energy in both spatial and loss limits. 

We can attribute the production of UHECRs in LLAGN based on the fact that in LLAGN the 
jet power can be magnetically provided by the BH spin-down power, which is dependent on the 
square of the strength of the BH magnetic field (Pj ~ B^). Thus, the particle acceleration regions 
can be sustained by the magnetic energy extraction from spinning BHs, where the strength of 
the magnetic field at the sites of particle acceleration scales with the maximum value of the BH 
magnetic field, which is ~ 10"^ gauss for a BH of 10^ M©, and where the maximum particle energy 
is given by Eqs. I47]andl49l 

The model proposed here is relevant for many nearby radio galaxies. Since their accretion 
rate is very low, the BH spin-down power, which is a better option than the accretion power, can 
sustain the required conditions for particle acceleration along the jets. Our predictions can be used 
in various propagation codes for UHECRs, which typically use the maximal energy and flux of 
the cosmic rays scaled relative to a canonical radio galaxy (e.g., columns 8 and 11 in Table 2), 
as well as the slope of the CR spectrum. The only difficulty is that the predicted quantities scale 
a slightly differently if some of the UHECRs are heavy nuclei at the origin. That is because for 
heavy nuclei, the energy scales up with the atomic number. For a middling FIR/Radio ratio ('FIR' 
is the flux density at 60 /i and 'Radio' is the flux density at 5 GHz), as Cen A has, a starburst may 
have provided all the necessary pre-accelerated nuclei. Nonetheless, this is not the case for M87. 

More information on radio data and BH mass from other nearby AGN are required to develop 
an improved statistical model which may answer the question of whether the flux of the UHECRs 
is produced by many weak radio galaxies or a few strong radio galaxies. 
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Table 1 : The jet parameters for M87 and Cen A and the corresponding estimation of the maximum particle energy 
(spatial and loss limits), luminosity, and flux of the UHECRs, were we assume that the particles are protons. 



Source 


7} 




e 




M 


rSGHz 
core 


^max 


ploss 
max 




^CR 


^jet 






(o) 


(o) 


(Mpc) 


(xlO'^Mo) 


(mJy) 


(10^° eV) 


(102" eV) 


(erg s"') 


(erg s"'cm"2) 


(erg s-') 


M87 


6 


10 


19 


16 


3 


2875.1 


8.61 


2.50 


2.03 X 10^-^ 


7.03x10-'*' 


~ 10^5 


Cen A 


2 


65 


5 


3.5 


0.055 


6984 


1.16 


2.17 


3.16 X 10^2 


2.28x10-'^ 


~ 10« 



NOTE: Col. (1) Source name; Col. (2) Jet Lorentz factor; Col. (3) Angle to the line of sight; Col. (4) Jet semi-opening angle; Col. (5) Distance to the source; Col. (6) BH 
mass; Col. (7) Core flux density at 5 GHz; Col. (8) maximum particle energy (spatial limit); Col (9) maximum particle energy (loss limit); Col. (10) UHE CR luminosity for a 
po wer-law index of paiticle distri bution p = 2; Col. (11) UH ECR flux; Co l (12) Esti mated energy flow alon g the jet Ljet. REFERENCES: J et Lorentz factor: 'Bir etta et al.l l ll999h 
and lMeisenheimer & etaTI )2007l) ; Angle to the Une of sight: iBiretta et afl i 19991) and lTingav & et al.l il998|); Distance to th e so urce for M87 iMacri & e t al. ( 1999) and for Cen A 
we as sume a distance Ds = 3 .5 Mpc to be consistent with the BH mass determination; BH mass: 'Macchetto & et alj il997h and lCappeUari & et alj i2009i) ; Core flux density at 5 
GHz: iNagar et al.l feOOll) and lSlee & et"ai] (1994;); Estimated energy flow along the jet Lj^O _Whvsons & A ntonucc r )2003l) . The value of the quantities infeiTed from observational 
data are the median ones, which we use in our estimations. 



Table 2: UHECR predictions for a complete sample of 29 steep spectrum sources (|Biermann & et al.ll2008l: ICaramete 
201lh . 



D 

(Mpc) 

(2) 



M 

(xlO'^Mo) 
(3) 



rSGHz 
core 

(mJy) 
(4) 



(10^0 eV) 
(5) 



E 

(6) 



max/ 
sp,M87 
max 



?loss 



(lO^o eV) 
(V) 



Floss / 
max' 
r^loss,M87 
■^max 



(8) 



^CR 

(erg s~') 
(9) 



(erg s~^ cm" 
(10) 



Fcr/ 

Z7M87 
^CR 

(11) 



Source 



_(1) 

ARP 308 69J Ol 67^8 L57 OTS 045 018 8.39x10^^ hSlxlO^ OlF 

CGCG 114-025 67.4 0.19 443.39 2.16 0.25 0.63 0.25 4.70x10^4 9.15x10"'° 1.30 

ESO137-G006 76.2 0.92 631.32 4.76 055 1.38 0.55 3.19x10^4 4.87x10"!° 069 

IC4296 54.9 1 214 4.97 057 1.44 0.57 3.75x10^^ 1.10xlO"i° 015 

IC5063 44.9 0.2 29 2.22 025 0.64 0.25 6.49x10^^ 2.84x10"" 0.04 

NGC0193 55.5 0.2 40 2.22 0.25 0.64 0.25 1.58x10^^ 4.55x10"" 0.06 

NGC0383 65.8 0.55 92 3.68 0.42 1.07 0.42 3.19x10^^ 6.52x10"" 0.09 

NGC 1128 92.2 0.2 39 2.22 0.25 0.64 0.25 5.20x10^3 5.41x10"" 0.07 

NGC 1167 65.2 0.46 243 3.37 039 0.98 0.39 l.BxlO^^ 2.36x10"^° 033 

NGC 1316 22.6 0.92 26 4.76 055 1.38 0.55 3.76x10^' 6.51x10"'^ 0.01 

NGC 1399 15.9 0.3 10 2.72 031 0.79 0.31 1.12x10^1 3.94x10"^^ o.005 

NGC2663 32.5 0.61 160 3.88 0.45 1.13 0.45 1.06x10^3 8.89x10" 012 

NGC3801 50 0.22 635 2.33 0.27 067 0.27 3.19x10^^ 1.12x10""^ 1.60 

NGC3862 93.7 0.44 1674 3.29 0.38 0.96 0.38 2.83x10^^ 2.85x10""^ 4.06 

NGC4261 16.5 0.52 390 3.58 0.41 1.04 0.41 6.80x10^^ 2.20x10"^° 0.31 

NGC4374 16 1 168.7 4.97 0.57 1.44 0.57 1.46x10^^ 5.05x10" 0.07 

NGC4486 16 3 2875.1 8.61 1 2.50 1 2.03x104^ 7.03x10"^° 1 

NGC4651 18.3 0.04 15 0.99 Oil 0.28 Oil 1.05x10^^ 2.78x10"" 0.03 

NGC4696 44.4 0.3 55 2.72 0.31 0.79 0.31 1.02x10^^ 4.59x10"" 0.06 

NGC5090 50.4 0.74 268 4.27 0.49 1.24 0.49 4.94x10^^ 1.72x10"'° 0.24 

NGC5128 3.5 0.055 6984 1.16 013 0.33 013 2.53x10^^ 1.82x10"*^ 25.95 

NGC5532 104.8 1.08 77 5.16 0.6 1.50 06 4.92x10^^ 3.96x10"" 0.05 

NGC5793 508 014 95.38 1.86 021 0.54 0.21 4.67x10^^ 1.60x10"'° 0.22 

NGC7075 72.7 0.25 20 2.48 0.28 0.72 0.28 l.BxlO^^ 1.89x10"" 0.02 

UGC 01841 84.4 01 182 1.57 018 0.45 018 4.34x10^^ 5.39x10"'° 0.76 

UGC 02783 82.6 0.42 541 3.22 0.37 0.93 0.37 5.58x10^^ 7.23x10"'° 1.02 

UGC 11294 63.6 0.29 314 2.67 031 0.78 0.31 2.01x10^^ 4.39x10"'° 062 

VV201 66.2 01 88 1.57 018 0.45 018 1.01x10^^ 2.04x10"'° 029 

WEIN45 84.6 0.27 321.6 2.58 0.3 0.75 0.3 4.31x10^^ 5.33x10"'° 0.75 



NOTE: Col. (1) Source name; Col. (2) Distance to the source; Col. (3) BH mass; Col. (4) Core flux density at 5 GHz; Col (5) Maximum 
particle energy (spatial limit); Col. (6) Maximum particle energy relative to that of M87; Col. (7) UHECR luminosity; (8) UHECR flux; (8) UHECR 
flux relative to that of M87. The value of the quantities inferred from observational data are the median ones, which we use in our estimations. 
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